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Abstract 

A search for antiproton decay has been made at the Fermilab Antiproton 
Accumulator. Limits are placed on six antiproton decay modes which contain 
a final-state muon. At the 90% C.L. we find that Tp/B(p — > fi~~ 7) > 5.0 x 
10 4 yr, T ¥ /B(p -> /i-vr ) > 4.8 x 10 4 yr, T f /B(p -> ^rj) > 7.9 x 10 3 yr, 
T ¥ /B(p -» ^-77) > 2.3 x 10 4 yr, T ¥ /B(p -> fi~K° s ) > 4.3 x 10 3 yr, and 
T f /B(p u~Kl) > 6.5 x 10 3 yr. 



PACS Numbers: 13.30Ce, 11.30.Er, 11.30.Fs, 14.20Dh 



1 Submitted to Phys. Rev. Lett. 



The CPT theorem requires that the proton and antiproton (p) lifetimes are 
equal. Searches for proton decay have yielded lower limits on the proton lifetime 
r p > O(10 32 ) yr 0. A search for p decay with a short lifetime (r^ << r p ) tests both 
the CPT theorem and the intrinsic stability of antimatter. 

Currently, the most stringent lower limits on Tp- can, in principle, be obtained 
from a comparison of recent measurements of the cosmic ray antiproton flux with 
predictions based on expectations for secondary production of antiprotons in the 
interstellar medium. The agreement between the observed and predicted rates implies 
that Tp is not small compared to T/7, where T is the p confinement time within the 
galaxy (~ 10 7 yr) and 7 is the Lorentz factor for the observed antiprotons. However, 
to extract a limit on Tp requires careful consideration of the relationship between the 
interstellar p flux and the flux observed at the Earth, and an adequate estimate of 
the systematic uncertainties on the predictions. When completed these calculations 
are expected to yield a limit in the range Tp > 10 5 -10 6 yr 0]. This indirect limit 
would not be valid if current models of p production, propagation, and interaction in 
the interstellar medium are seriously flawed by, for example, our lack of knowledge of 
the nature of most of the matter in the galaxy. 

Laboratory searches for p decay have so far provided less stringent limits on 
Tp. However, these limits do not suffer from large model dependent uncertainties. 
The most stringent published laboratory limits on p decay have been obtained from 
(a) measurements of the containment lifetime of ~ 1000 antiprotons stored in an ion 
trap, yielding |§ Tp > 3.4 months, and (b) prior searches at the Fermilab antiproton 
accumulator for explicit p decay modes with an electron in the final state, yielding || 
T ¥ /B(p — > e~7) > 1848 yr, Tp/B(p — > e~7r°) > 554 yr, T ¥ /B(p — > e^ri) > 171 yr, 
Tp/B(p -> e~K s ) > 29 yr, and T ¥ /B{p e~K° L ) > 9 yr. 

Angular momentum conservation requires that a decaying p would produce a 
fermion (electron, muon, or neutrino) in the final state. In this paper we report 
results from the first search for antiproton decay modes with a muon in the final 
state. The results have been obtained using the APEX detector at the Fermilab 
Antiproton Accumulator. The experiment was designed to search for decay modes 
with an electron or muon in the final state. The results from the search for p —> e~+X 
will be described elsewhere ||. 

A full description of the APEX detector can be found in Ref. 0. In the following 
we give a brief description of the main detector components relevant to the analysis 
described in this paper. A schematic of the APEX detector is shown in Fig. [I]. The 
detector, which was located in a straight section of the Accumulator ring, consisted 
of: (i) A 3.7 m long evacuated decay tank operated at 10~ n Torr. The downstream 
section of the tank consists of a 96 cm diameter cylinder with a 1.2 mm thick stainless 
steel window, (ii) A movable tungsten wire target at the upstream end of the tank 
(z = 0, where in the APEX co-ordinate system z is measured in the direction of the 
antiproton beam). The target could be inserted into the beam halo to produce a 
localized source of particles for aligning and calibrating the detector, (iii) Horizontal 
and vertical scintillation counters arranged around the 10 cm diameter beam pipe 
immediately upstream of the tank and target. The counters covered a 1 x 1 m 2 
area normal to the beam direction, and were used to veto tracks from upstream 



2 



interactions, (iv) Three planes of horizontal and three planes of vertical scintillation 
counters downstream of the tank. The last planes of horizontal and vertical counters 
were downstream of a 2.3 radiation length lead wall, providing a preradiator to aid 
in identifying electrons and photons. The remaining counter planes were upstream 
of the lead, and provided pulse height information used to count the number of 
charged particles in an event (dE/dx counters), (v) Three planes of horizontal and 
three planes of vertical 2 mm diameter scintillating fibers downstream of the tank 
and upstream of the preradiator lead. These detectors were used for tracking, and 
provided three space points along the track trajectory with typical residuals of 620 jum 
in the directions transverse to the beam direction. This enabled the origin of tracks 
emerging from the decay tank to be located with a precision a z ~ 13 cm along 
the beam orbit. The measured average single-hit efficiency for the tracker planes is 
(89.2 ± 1.5)% for tracks in the event samples described in this paper. The measured 
track reconstruction efficiency, based on a sample of events that are consistent with 
having one minimum-ionizing particle passing through the dE/dx and muon telescope 
(see below) counters, is (90 ± 7)%. (vi) A lead-scintillator sampling electromagnetic 
calorimeter [7| constructed from 144 rectangular 10 x 10 cm 2 modules that are 14.7 
radiation lengths deep. The modules were arranged in a 13 x 13 array with 6 modules 
removed from each of the four corners, and the central module removed to allow 
passage of the beam pipe. The calorimeter was calibrated by measuring the response 
to minimum ionizing tracks and reconstructing tt° — > 77 and 77 — > 77 mass peaks 
using data samples recorded with the calibration target inserted in the beam halo. 
The measured mass peaks had fractional r.m.s. widths given by cr m /m ~ 0.25 for 
cluster pairs in the energy range of interest for the experiment. This mass resolution 
is dominated by the energy resolution of the calorimeter, (vii) A tail catcher (TC) 
downstream of the calorimeter consisting of a 20 cm deep lead wall followed by planes 
of horizontal and vertical scintillation counters, (viii) A limited-acceptance muon 
telescope (MT), 10 nuclear interaction lengths deep, located downstream of the TC, 
and aligned to point towards the center of the decay tank. The MT consists of a 
sandwich of five iron plates and five 30 x 30 cm 2 scintillation counters. 

The APEX experiment took data when there were typically 10 12 antiprotons 
circulating in the 474 m circumference accumulator ring operating with a central p 
momentum of 8.90 ± 0.01 GeV/c (7 = 9.54 ± 0.01). A measure of the sensitivity of 
the APEX data sample is given by: 

S = - / Np(t)dt = (3.31 ±0.03) x 10 V, 

7 J 

where Np (t) is the number of circulating antiprotons at time t, the integral is over 
the live-time of the experiment, and the uncertainty arises from the precision with 
which the time dependent beam current was recorded. 

To search for muonic p decays occurring as the beam particles traverse the decay 
tank, data were recorded with a muon trigger that required a coincidence between 
at least two of the five MT scintillation counters. These triggers were eliminated if 
they were in coincidence with a signal in one or more of the upstream veto counters 
indicating the presence of an interaction upstream of the decay tank. This loose 
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trigger resulted in 1.2 xlO 6 events being recorded. These events predominantly arise 
from interactions of the p beam with the residual gas in the decay tank or with 
material surrounding the beam. The coincident MT counter signals are then produced 
by traversing muons coming from charged pion decays, and by hadronic showers not 
contained in the calorimeter and TC. 

Offline, after final calibration of the scintillation counters, the upstream veto 
counter requirement was re-applied using a more stringent threshold. This reduced 
the data sample to 1.1 xlO 6 events. Further requirements were then imposed to select 
events containing a candidate energetic muon that traverses the MT and originates 
from within the decay tank. These requirements were (i) that at least four of the 
five MT counters be above threshold (4.2 x 10 4 events), and (ii) the presence of one 
and only one scintillating fiber track that extrapolates to the MT counters within the 
expected uncertainty due to multiple scattering, and also extrapolates back to the 
beam orbit with a point of closest approach within the fiducial volume of the decay 
tank (0 < z < 375 cm) and with an impact parameter less than 1 cm (416 events). 

To further suppress backgrounds additional requirements can be imposed on the 
event topology and kinematics. We begin by considering the process p — > /i~7, which 
would result in events in which an energetic photon is produced coplanar with the 
muon (i.e. traveling within the plane defined by the muon and the incoming beam 
direction). After requiring at least one neutral cluster || in the calorimeter (209 
events) that is coplanar with the muon candidate (±5°), the data sample is reduced 
to 14 events. The observed neutral cluster energy distribution for these events is 
compared in Fig. Q with the predicted distribution obtained using the GEANT || 
simulation described below, and corresponding to p — > /i~7 decay with a lifetime 
Tp/B(p — > fi~j) = 5000 yr. The observed distribution peaks at low cluster energies, 
with a tail extending to approximately 2 GeV. In contrast, the predicted distribution 
for p — > /i~7 decays peaks at about 3.5 GeV, with only 4.4% of the simulated events 
having cluster energies less than 2 GeV. We conclude that there is no evidence for 
a signal. To eliminate the background we therefore require that the neutral cluster 
energy exceed a minimum value E min , and choose E min = 3 GeV. No candidate events 
remain. We note that although the choice of E m i n is somewhat arbitrary, the final 
limit that we obtain on the decay p — > pT^ is not very sensitive to small changes in 
Emi n . The resulting limit on Tp/B(p —> pT'y) is given in years by 

r ¥ /B{p - yrj) > - - 1 — / N v (t) dt = (3.31 ± 0.03) x 10 9 , (1) 

7 -'"max l^max 

where N max = 2.3 is the 90% C.L. upper limit on the observation of N = events, 
and e is the fraction of decays taking place uniformly around the accumulator ring 
that would pass the trigger and event selection requirements. To take account of a e , 
the systematic uncertainty on e, we use the prescription given in Ref. 0, giving at 
90% C.L. 

N max = 2.3 x (1 + 2.3^/2), (2) 

where a r = cr e /e. 

The GEANT Monte Carlo program has been used to simulate the detector re- 
sponse and calculate e. The detector simulation includes a full description of the 
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detector geometry, and correctly describes the calorimeter, tracker, and relevant scin- 
tillation counter responses (dE/dx, MT) measured using calibration data, together 
with the measured performance of the muon trigger. Further details can be found in 
Ref. ||. The efficiency e was obtained by generating 10 5 p — > /i~7 decays uniformly 
along the beam orbit within the decay tank. To a good approximation the geometri- 
cal efficiency of the detector and trigger is negligible for decays occurring outside of 
the tank. We obtain e = (3.7 ± 0.9) x 10~ 5 , where the uncertainty on e arises from 
the systematic uncertainties on the trigger calibration, calorimeter energy scale, and 
track efficiency. The trigger and calorimeter scale uncertainties yield contributions to 
a e of ±18% and ±15%, respectively, and have been evaluated by analyzing GEANT 
Monte Carlo p —> /i~7 samples in which the simulated trigger and calorimeter scales 
have been changed by ±1 a. The overall systematic uncertainty on e has been calcu- 
lated by adding these contributions in quadrature with the uncertainty on the track 
efficiency (±7%). Inserting the calculated e and a e into Eqs. 1 and 2 we obtain the 
result T ¥ /B{p -» /jl'j) > 5.0 x 10 4 yr (90% C.L.). 

Now consider the two-body decays p — > /i~7r° and p — > fi~r] (with 7r°, rj — > 77), 
and the three-body decay p — > /i~77- These decays would result in events with 
one or two neutral clusters observed in the calorimeter, where the one-cluster events 
occur when the two photon-showers are not spatially resolved in the calorimeter or 
when one of the photons is outside of the calorimeter acceptance. In addition, the 
77 may also decay into more complicated final states producing further clusters in 
the calorimeter. To optimize the search for fi~ir and fi~r] final states we divide the 
209 events described previously that have a muon candidate plus one or more neutral 
clusters into two subsamples, namely a one-cluster sample containing 104 events, and 
a multi-cluster sample containing 105 events. The clusters in the one-cluster sample 
are required to be coplanar with the muon (±10°), and the cluster-pairs formed by 
the two highest energy clusters in the multi-cluster events are also required to be 
coplanar with the muon (±10°). These requirements reduce the samples to 13 one- 
cluster and 10 multi-cluster events. We next require that the multi-cluster events, 
which would contain 2 photons if they were genuine p — > ^~tt°, /i - ^, or /i~77 decays, 
have a preradiator signal above a threshold of 0.5 x minimum-ionizing pulse height. 
Only 3 multi-cluster events satisfy this requirement. Hence we are left with 16 events 
with > 1 cluster events (/i ± > 1 cluster) for further analysis. 

Under the hypothesis that the observed (/ry) and (^77) systems arise from the 
decay of a beam particle, the mass of the parent particle can be computed from the 
measured muon direction and the directions and energies of the neutral clusters, using 
the constraint that the vector sum of the momentum components of the daughter 
particles transverse to the beam direction is zero. The resulting mass distribution 
for the remaining 16 (/x ± > 1 cluster) events is compared in Fig. |^ with predictions 
from the GEANT simulation for the decays (a) p — > /U~7r°, (b) p — > fi~r], and (c) 
p — > /U~77. The observed mass distribution peaks at low masses with a tail extending 
to approximately 0.7 GeV/c 2 . In contrast to this, the simulated signal distributions 
peak at the p mass, with only 17% (23%) [13%] of the simulated /i~7r° (/i~r/) [/i~77] 
decays resulting in reconstructed masses less than 0.7 GeV/c 2 . We conclude that there 
is no evidence for a signal. To eliminate the background we therefore require that the 
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reconstructed mass exceed a minimum value m min , and choose m min = 0.75 GeV/c 2 . 
No candidate events remain. We note that although the choice of m min is somewhat 
arbitrary, the final limits that we obtain on the decays p — *> /i~vr°, fi~f], and /i~77 are 
not very sensitive to small changes in m m i n . The calculated overall efficiencies e for 
these decays to pass our trigger and analysis requirements are (3.6±0.9) x 10 -5 for the 
/i~7T° mode, (6.1 ± 1.9) x 10~ 6 for the \x~r\ mode, and (1.8 ± 0.5) x 10 -5 for the /i~77 
mode. Substituting the values into Eqs. 1 and 2 yields the limits Tp/B(p — > > 
4.8 x 10 4 yr, T ¥ /B(p -> ynj) > 7.9 x 10 3 yr at 90% C.L., and T ¥ /B(p -» ^-77) > 
2.3 x 10 4 yr at 90% C.L. 

Now consider the other possible two-body muonic decay modes of the antipro- 
ton, namely decays into the final states n~Kg, fJ,~K®, p,~p°, and pTuj. Detailed 
GEANT simulations have been made for these decay modes. The muons from the 
p~p° and f/Tu decay modes are predicted to have energies that are too low to pen- 
etrate the muon telescope. We therefore restrict ourselves to the fJTKg and \i~K\ 
modes. The calculated efficiencies e for these decays to satisfy the trigger and ei- 
ther the /i~7 or /i~7r° search criteria described previously are (3.3 ±1.0) x 10~ 6 for 
the fi~K§ mode, and (5.0 ± 1.5) x 10~ 6 for the p~K\ mode. Substituting these 
values into Eqs. 1 and 2 yields the limits Tp/B(p — > n~Kg) > 4.3 x 10 3 yr and 
T ¥ /B(p -> n~Kl) > 6.5 x 10 3 yr at 90% C.L. 

Finally, our results for the five muonic decay modes presented in this paper are 
summarized in Table 1, and a more comprehensive description of the analysis can be 
found in Ref. 0. 

The APEX experiment was performed at the Fermi National Accelerator Lab- 
oratory, which is operated by Universities Research Association, under contract DE- 
AC02-76CH03000 with the U.S. Department of Energy. 
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Table 1: Summary of lifetime limits for p decay to six muonic final states with calcu- 
lated efficiencies e for each mode. 



Decay Mode e r/B Limit (yr) 

(90 % C.L.) 

fi- +7 (3.7 ±0.9) x 1CT 5 > 5.0 x 10 4 

/i- + 7r° (3.6 ±0.9) x 10" 5 >4.8xl0 4 

H' + rj (6.1 ± 1.9) x 10- 6 > 7.9 x 10 3 

At" ±77 (1.8 ±0.5) x 10" 5 >2.3xl0 4 

At"±^ (3.3 ± 1.0) x 10" 6 >4.3xl0 3 

fT + Kl (5.0 ± 1.5) x 10~ 6 >6.5xl0 3 
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Figure 1: Schematic of the APEX detector. 
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Figure 2: Distribution of neutral cluster energies for the 14 events that pass the 
p — > /i~7 selection criteria described in the text (open histogram) compared with 
the predicted distribution for a signal corresponding to Tp/B = 5000 yr (hatched 
histogram). Limits are based on the region E > 3 GeV. 
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Figure 3: Distribution of event masses for the 16 events that pass the selection criteria 
described in the text (open histograms) compared with predictions from a GEANT 
simulation (hatched histograms) for the decays (a) p — > /i~vr°, (b) p — > 77, and 
(c) p — * // _ 77- The two entries in the lowest bin are events with calorimeter energy 
less than the 7r° rest energy. The predicted signal distributions are normalized to 
correspond to Tp/B = 5000 yr for the fi~n° mode, 1000 yr for the /J,~r] mode, and 
2000 yr for the /i~77 mode. Limits are based on the region Mass > 0.75 GeV. 

11 



